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I.  INTRODUCTION 


This  is  the  final  report  on  this  program  which  was  a  collaborative  effort  among  AeroChem 
(Principal  Investigator  H.  F.  Calcote),  Penn  State  (Principal  Investigator  M.  Frenklach)  and  Iowa 
State  (Principal  Investigator  R.  C.  Brown).  The  ultimate  objective  has  been  to  develop  a  quantitative 
model  of  soot  formation  in  flames  that  is  consistent  with  experimental  data.  The  specific  objectives 
of  this  collaborative  three-year  study  were:  (1)  to  delineate  the  relative  importance  of  the  neutral  free 
radical  and  ionic  mechanisms  of  soot  formation  in  flames;  (2)  to  determine  the  optimum  model  of  the 
total  soot  formation  process  based  upon  current  knowledge;  and  (3)  to  recommend  how  to  improve 
the  model  and  the  additional  experiments  necessary  to  clarify  any  discrepancies. 

Professor  Robert  Brown,  Iowa  State  University,  was  added  to  the  program  in  the  last  year  as 
a  collaborator  via  a  subcontract  from  AeroChem,  with  funds  which  were  added  to  this  program  for 
that  purpose.  He  brought  to  the  program  previous  experience  in  developing  a  code  for  modeling  ionic 
reactions  in  flames1,2  involving  the  solution  of  the  stiff  equations  involved  in  such  reactions.  Brown’s 
code  was  originally  developed  to  include  ions  as  well  as  neutrals.  It  handles  non-Arrhenius 
temperature  coefficients  in  the  ion-molecule  rate  coefficients;  in  the  code  used  by  Frenklach  this  has 
been  a  problem.  The  diffusion  coefficient  for  ions  and  electrons  is  included  via  Poisson’s  equation; 
in  the  code  used  by  Frenklach  this  is  handled  in  an  indirect  way.  Brown’s  code  also  permits  running 
the  neutral  species  separately  and  then  using  these  results  as  input  to  the  ion  code.  This  allows  the 
use  of  experimental  species  profiles  as  input,  so  that  only  the  part  of  the  mechanism  of  interest  can 
be  studied.  As  we  will  see  this  makes  a  major  difference.  The  validity  of  the  test  of  the  ionic 
mechanism  is  no  longer  dependent  upon  the  free  radical  mechanism  of  small  species.  We  have  already 
spent  considerable  effort  in  developing  this  part  of  the  mechanism  and  have  not  yet  had  a  good  test 
of  it.  As  we  will  see,  the  free  radical  mechanism  is  not  consistent  with  the  observed  small  species  in 
the  flame.  A  more  complete  set  of  experimental  neutral  species  profiles  was  made  available  to  us  by 
Vovelle.3 

The  AeroChem  and  Iowa  State  effort  are  presented  in  this  report;  the  Penn  State  effort  is 
covered  in  their  companion  report.4  AeroChem  is  responsible  for  development  of  the 
thermodynamics,  diffusion  coefficients,  reaction  mechanism,  and  reaction  rate  coefficients  for  the 
ionic  mechanism.  There  were  no  previous  sets  of  data  on  the  ionic  mechanism  to  draw  upon  when 
we  initiated  this  program,  two  papers  were  subsequently  published  by  Brown  and  Eraslan,1'2  in  which 
they  calculated  the  concentrations  of  ions  in  both  stoichiometric  and  fuel  rich  flames;  their 
calculations  compared  favorably  with  experiments.  For  the  sooting  flames  we  model  in  this  program, 
it  has  been  necessary  not  only  to  construct  the  ionic  mechanism  and  identify  the  rate  coefficients  but 
also  to  develop  neutral  mechanisms  and  to  organize  data  for  odd  number  neutral  carbon  compounds. 
Such  species  were  not  generally  utilized  in  the  free  radical  mechanisms  of  soot  formation.  Odd  carbon 
neutral  species,  however,  play  a  significant  role  in  the  ionic  mechanism  and  are  observed  in  relatively 
large  concentrations  in  sooting  flames. 

The  computer  modeling  results  are  compared  with  the  well-documented  acetylene/oxygen 
flame  (the  "standard  flame")  burning  on  a  flat  flame  burner  at  a  pressure  of  2.67  kPa  (20  Torr)  and 
a  linear  flow  rate  of  50  cm/s  in  the  unburned  gases.  In  our  previous  AFOSR  contract  work,5-7  we 
duplicated  the  Bittner  and  Howard8  burner  used  to  obtain  neutral  species  concentrations  and  we 
measured  ion  concentration  profiles  with  this  burner  so  that  the  ion  profiles  and  the  neutral  profiles 
would  be  from  the  same  system.  We  also  compared  the  data  obtained  by  a  number  of  other 
researchers  on  very  nearly  the  same  flame,  we  previously  presented  the  results  of  this  comparison5 
which  showed  amazing  agreement  among  several  laboratories.  This  data  base  forms  an  excellent 
experimental  standard  with  which  to  compare  the  computer  modeling  results. 
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Preliminary  computer  experiments  were  carried  out  for  a  shock  tube  to  obtain  some  feel  of 
how  the  ionic  mechanism  performed  in  the  computer  program,  and  to  work  out  some  of  the  details 
of  the  mechanism.  This  is  not  a  good  system  in  which  to  compare  the  two  mechanisms  because  there 
are  no  ion  profile  data.  Further,  at  the  time  of  this  preliminary  work  the  ionic  mechanism  was  in  a 
very  crude  state  of  development  so  the  results  are  of  little  value  in  so  far  as  comparing  the  ionic  and 
free  radical  mechanisms.  Frenklach  has,  however,  chosen4  to  interpret  these  results  as  indicating  "that 
the  ionic  mechanism  produces  polycyclic  aromatic  hydrocarbons  (PCAH)  at  a  significantly  lower  rate 
than  does  the  mechanism  composed  of  only  reactions  of  neutral  species."  Not  only  is  this  a  premature 
conclusion,  but  more  seriously,  Frenklach  continues  to  misinterpret  the  ionic  mechanism,  in  spite  of 
many  personal  explanations. 

The  basic  premise  of  the  ionic  mechanism  is  that  ions  grow  to  a  very  large  size  to  become 
incipient  soot  particles  or  to  be  neutralized  and  produce  large  neutrals,  500  to  1,000  u,  which  grow 
to  become  incipient  soot  particles.  Frenklach  interprets  the  ionic  mechanism  as  producing  small 
PCAH  neutrals  and  then  compares  the  concentration  of  such  species  as  benzene  and  acenaphthalene 
produced  by  the  free  radical  and  ionic  mechanisms.  Frenklach  states  in  his  abstract,4  "the  formation 
of  polycyclic  aromatic  hydrocarbons,  the  precursors  to  soot,  via  the  ionic  reaction  pathway  is  slower 
than  via  the  pathway  involving  neutral  species";  neutral  polycyclic  aromatic  hydrocarbons  do  not  play 
a  role  (except  as  by-products)  in  the  ionic  mechanism.  The  production  of  such  small  neutrals  is  a 
result  of  ion-molecule  reactions  which  remove  ions  from  the  main  stream,  and  are  thus  unwanted  by¬ 
products.  The  basic  question  which  should  be  addressed  is  what  is  the  flux  of  carbon  (atoms)  through 
the  ionic  reaction  scheme  vs.  through  the  radical  reaction  scheme  of  leading  to  soot;  neutrals  for  the 
free  radical  mechanism  and  ions  for  the  ionic  mechanism. 

At  the  present  state  of  this  research,  it  would  be  premature  to  decide  between  the  two 
mechanisms.  We  present  some  new  evidence  which  favors  the  ionic  mechanism  over  the  free  radical 
mechanism.  Reference  9  summarizes  some  of  the  earlier  evidence  for  the  ionic  mechanism.  There 
is  no  equivalent  evidence  for  the  free  radical  mechanism.  Consistent  with  the  contract  objectives,  we 
identify  work  which  needs  to  be  done  to  understand  the  relative  role  of  ions  and  neutrals  in  soot 
inception. 


II.  STATEMENT  OF  WORK 

A.  Organize  relevant  data  on  ion-molecule  reactions,  thermochemistry,  electron  attachment,  ion 
recombination,  and  ion  and  electron  diffusion  to  be  used  in  the  computer  code  simulation  by 
Penn  State. 

B.  Determine  the  hydrogen  atom  concentration  in  the  "well-studied"  C2H2/02  flame  and 
determine  if  this  concentration  exceeds  the  thermal  equilibrium  concentration. 

C.  Analyze  the  computer  simulation  data  obtained  by  Penn  State  and  compare  the  results  with 
available  experimental  data  to  determine  the  major  chemical  pathways  to  incipient  soot  and 
to  simplify  the  computer  model. 

D.  Organize  relevant  data  on  the  elementary  steps  involved  in  the  growth  of  incipient  soot  to  soot 
particles,  including  growth  by  molecular  addition,  coagulation,  and  oxidation. 

E.  Analyze  the  computer  simulation  data  obtained  by  Penn  State  using  the  extended  model  and 
compare  the  results  with  available  experimental  data  to  determine  the  major  pathways  to 
particulate  soot  and  to  determine  how  to  alter  the  model  to  make  it  more  in  conformity  with 
experimental  results. 
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F.  Review  the  literature  and  choose  flame  experiments  with  which  to  compare  the  model 
developed  above. 

G.  Analyze  the  computer  simulation  data  obtained  by  Penn  State  and  compare  the  results  with 
the  experimental  data  to  determine  the  general  applicability  of  the  model  and  to  recommend 
what  is  required  to  improve  the  model  and  what  additional  experiments  are  necessary  to 
clarify  any  discrepancies. 

The  Statement  of  Work  for  the  subcontract  with  Brown  at  Iowa  State  was: 

A.  Incorporate  reaction  mechanism  and  thermochemical  data  developed  at  AeroChem  Research 
Laboratories  in  the  TRANSEQI  computer  code. 

B.  Assess  the  utility  of  determining  ambipolar  diffusion  by  numerical  solution  of  the  Poisson 
equation  using  a  sample  ion  model. 

C.  Perform  baseline  simulations  of  a  test  flame. 

D.  Perform  numerical  simulations  of  the  standard  acetylene/oxygen  flame. 


III.  THERMODYNAMICS 


The  required  thermodynamic  data  for  most  of  the  ions  and  some  of  the  neutral  species  were 
either  collected  from  the  literature  or  calculated.  Available  compilations  of  data  for  ions,  e.g.,  Refs. 
10-12  contain  very  few  odd  number  carbon  atom  ions;  experimentally  these  species  dominate  in  the 
flame  mass  spectra.  The  compilations  also  do  not  include  isomers,  which  are  of  interest. 
Thermodynamic  data  for  these  species  were  generated  as  described  below. 

Data  on  a  number  of  species,  e.g.,  vinyl  radical,  C3H2  isomers,  C2H3+  isomers,  C3H3+  isomers, 
c7h5+,  c8h7+,  c10h9+,  C13H8+,  and  several  larger  ions,  were  recalculated,  based  upon  new 
thermodynamic  information.  New  calculations  were  made  for  94  ions  and  their  isomers.  The  results 
for  many  of  these  were  presented  in  Ref.  13  and  will  not  be  included  here. 

The  thermodynamic  quantities,  Cp°  and  S°  for  hydrocarbons  and  hydrocarbon  ions,  were 
calculated,  when  sufficient  information  was  available  from  vibrational  frequencies  and  moments  of 
inertia  using  statistical  mechanical  methods,  or  when  not  available,  following  Benson’s 
thermochemical  methods  of  group  additives.  Enthalpies  of  formation  at  298  K,  AHf°  298,  were 
obtained  either  directly  from  the  literature  or  calculated  from  experimental  proton  affinities  and  the 
corresponding  neutral  molecule’s  enthalpy  of  formation. 

Unfortunately,  there  are  several  neutral  molecules  for  which  the  thermodynamic  data  are 
uncertain.  This  is,  for  example,  the  case  with  diacetylene,  an  important  reactant  in  the  ionic 
mechanism.  The  available  heats  of  formation  at  298  K  range  from  439  to  473  kJ/mol. 14-18  We  have 
used  440  KJ/mol14  because  the  data  in  this  reference  were  used  previously  in  several  thermochemical 
calculations  with  large  ions.  This  uncertainty  will  be  manifest  in  the  flame  simulations  when  rates 
for  reverse  chemical  reactions  are  calculated  from  thermodynamics  and  forward  reaction  rate 
coefficients.  For  a  few  cases,  such  as  for  the  cyclopropenyl  ion,  C3H3+,  another  important  participant 
in  the  ionic  mechanism,  there  was  sufficient  experimental  information  to  employ  the  more  accurate 
statistical  mechanical  methods  to  calculate  Cp°  and  S°. 
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The  inaccuracy  of  thermodynamic  data  continues  to  contribute  a  major  question  in  any 
reaction  kinetics  scheme. 


IV.  REACTION  MECHANISM  AND  REACTION  COEFFICIENTS 


The  overall  mechanism  is  shown  in  Fig.  1.  The  ionic  path  starts  when  electronically  excited 
CH  (CH*)  reacts  with  an  oxygen  atom  to  produce  HCO+,  which,  through  a  series  of  ion-molecule 
reactions,  produces  C3H3+,  an  ion  which  is  observed  in  large  concentrations  in  fuel  rich  and  sooting 
flames.  This  ion  then  reacts  with  any  of  six  small  neutral  species,  indicated  in  Fig.  1 ,  to  produce 
larger  ions  which  continue  to  grow  through  a  series  of  ion-molecule  reactions  with  the  same  six 
species  producing  larger  and  larger  ions.  To  simplify  the  number  of  reactions  that  have  to  be  handled 
we  have  only  considered  C2H2,  C4H2  and  C3H4  in  the  reaction  set.  Simultaneously  as  the  ions  grow, 
they  are  neutralized  at  a  rate  which  increases  with  increasing  mass,  producing  "neutral  by-products," 
polycyclic  aromatic  hydrocarbons,  that,  of  course,  can  continue  to  grow  to  soot  through  the  "free 
radical  mechanism."  These  neutral  reactions  have  not  yet  been  incorporated  into  the  mechanism;  they 
probably  play  only  a  small  role  except  for  very  large  molecules. 

The  complete  reaction  mechanism  with  reaction  rate  coefficients  and  free  energies  of  reaction 
at  1750  K  is  presented  in  Table  I  as  it  presently  exists.  This  set  of  128  reactions  has  been  chosen  from 
about  400  reactions. 

In  developing  this  mechanism  only  ionic  species  which  have  been  observed  in  sooting  flames 
have  been  included  and  all  ions  observed  in  the  flame  have  been  accounted  for.  This  is  a  more 
stringent  constraint  than  has  ben  applied  to  the  free  radical  mechanism. 


A.  EXCITED  STATE/CHEMIIONIZATION 

•  *  •  •  •  • 

Specific  reactions  which  produce  excited  CH,  (CH  ),  important  in  the  chemiionization  process, 

and  the  subsequent  chemiionization  reactions  are  presented  in  Table  I.  A.  Probably  the  two  most 

important  reactions  in  this  set  (a  sensitivity  analysis  has  yet  to  be  analyzed),  are  Reactions  (3)  and  (8). 

The  rate  coefficient  for  Reaction  (3)  was  estimated  by  W.  Gardiner19  and  the  rate  coefficient  for 

Reaction  (8)  was  measured  by  Cool  relative  to  the  rate  coefficient  for  Reaction  (7).20  The  reaction 

coefficient  for  Reaction  (7)  is  from  a  bimolecular  quantum  RRK  theory  calculation  by  Westmoreland 

which  is  in  agreement  with  experiment.3,14 


B.  ION-MOLECULE  REACTIONS 


When  available,  experimental  rate  coefficients  were  used,  but  these  are  available  only  for 
small  ions.  In  general,  experimental  rates  are  very  close  to  the  rate  calculated  by  the  average  dipole 
orientation,  ADO,  theory16 


k 


27re 

T1* 


+  C/rD 


2  ] 
7rkT 


(1) 


where,  n  is  the  reduced  mass,  a  is  the  polarizability  of  the  neutral  reactant,  C  is  a  locking  constant 
determined  from  experimental  data,  and  /iD  is  the  dipole  moment  of  the  neutral  reactant.  For 
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TABLE  I 

REACTION  MECHANISM  AND  REACTION  RATE  COEFFICIENTS 
FOR  THE  IONIC  MECHANISM  OF  SOOT  FORMATION 


k  =  A  Tn  e‘£/RT  (K,  kJ,  moles,  cm2,  s) 


A.  EXCITED  STATE/CHEMI IONIZATION  REACTIONS 


No. 

R 

E  A  C  T  I  0 

N 

.  A 

n 

E 

* 

1 

C2 

+ 

OH 

CH 

* 

+ 

CO 

3.4E+12 

0 

0 

2 

3: 

CM 

O 

+ 

0 

CH 

* 

+ 

CO 

7.1E+11 

0 

0 

3 

C,H 

+ 

°2 

CH 

+ 

C02 

4.5E+15 

0 

105 

4 

CH 

+ 

M 

-♦ 

CH 

+ 

H 

4.0E+10 

0.5 

0 

5 

* 

CH 

* 

+ 

°2 

-♦ 

CH 

+ 

°2 

2.4E+12 

0.5 

0 

6 

CH 

CH 

1.7E+06 

0 

0 

7 

CH 

+ 

0 

-♦ 

HCO+ 

+ 

e 

1.4E+10 

0 

2 

8 

* 

CH 

+ 

0 

- 

HCO+ 

+ 

e 

4.8E+14 

0 

0 

B. 

No. 

ION-MOLECULE 

R  E  A  C  T  I  0 

REACTIONS 

N 

a 

AG( 1750) 

500 

Forward  Rate  Coefficient,  k 
cm  /mole/s 

1000  1500  2000 

1 

ch3+  + 

c2h2 

= 

H3C3  + 

H2 

-213.1 

8.3E+14 

2 

ch3+  + 

c4h2 

= 

¥s+ 

-290.9 

1.0E+15 

3 

h3°+  + 

C3H2 

= 

H3C3+  + 

h2o 

-253.9 

8.5E+14 

4 

c2h3+  + 

C2^2 

= 

V5+ 

-96.7 

7.1E+14 

5 

C2H3+  + 

c4h2 

= 

C6«5+ 

-141.3 

8.2E+14 

6 

HC0+  + 

ch2 

= 

CH/  + 

CO 

-203.1 

5.2E+14 

7 

HC0+  + 

h2o 

= 

h30+  + 

CO 

-107.6 

1.9E+15 

8 

HC0+  + 

C2H2 

= 

C2H3+  + 

CO 

-75.0 

6.9E+14 

9 

HC0+  + 

CgHg 

= 

H3C3+  + 

OH 

32.2 

1.5E+13 

2.4E+13 

3.3E+13 

4.1E+13 

10 

HC0+  + 

C3H2 

= 

H3C3+  + 

CO 

-361.5 

7.5E+14 

11 

HC0+  + 

C3H4 

= 

H3C3+  + 

h2  +  CO 

-271.0 

7.9E+14 

12 

c3H3+  + 

M 

= 

H3C3+  + 

M 

-64.6 

5.0E+14 

13 

C3H3+  + 

C2H2 

- 

c5h2h+  + 

h2 

6.8 

6.5E+14 

6.5E+14 

6.5E+14 

4.5E+14 

14 

C3H3+  + 

c2h2 

= 

h5c5+ 

-135.6 

6.5E+14 

15 

C3H3+  + 

c4h2 

= 

C5H2H+  + 

¥2 

0.0 

7.4E+14 

7.4E+14 

7.4E+14 

5  1E+14 

16 

C3H3+  + 

c4h2 

= 

c7h4h+ 

-67.3 

7.4E+14 

17 

C4H5+  + 

CgH2 

= 

C6H5+  t 

H, 

-37.8 

6.2E+14 

18 

C5H2H+  + 

C2H2 

- 

c7h4ii+ 

-67.2 

6.0E+14 

19 

c5h2h+  + 

C3H4 

s 

C6H5+  + 

c2h2 

-178.9 

6.5E+14 
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No. 

REACT 

I  0 

N 

500 

1000 

1500 

2000 

20 

H5C5+  +  C2H2 

= 

H7C7+ 

-35.7 

6.0E+14 

21 

H5C5+  +  C4H2 

= 

¥o+ 

-30.4 

6.5E+14 

22 

C6H5+  +  CH2 

= 

c7h4h+  + 

H2 

-202.5 

4.7E+14 

23 

C6H5+  +  C2H2 

= 

-42.0 

5.8E+14 

24 

c7h/+  h2 

- 

»7C7+ 

-101.0 

8.9E+14 

25 

c?h4h+  +  C2H2 

= 

c9h7+ 

-17.6 

5.7E+14 

26 

c?h4h+  +  C2H2 

= 

h7c9+ 

-98.7 

5.7E+14 

27 

H?C7+  +  C2H2 

a 

h7c9+  + 

h2 

2.3 

5.7E+14 

5.7E+14 

5.7E+14 

4.6E+14 

28 

H7C7+  +  C2H2 

= 

C9H7+  + 

H2 

83.4 

7.7E+06 

5.2E+10 

1.2E+12 

5.8E+12 

29 

C8H7+  +  C2H2 

a 

C10H9+ 

-49.3 

5.6E+14 

30 

C8H7+  +  C2H2 

a 

H9C10+ 

81.0 

5.6E+14 

31 

C8H7+  +  C3H4 

a 

¥ll+  + 

H2 

-19.9 

6.0E+14 

32 

CgH7+  +  M 

= 

H7C9+  + 

H 

-81.1 

4.5E+14 

33 

C9H7+  +  C2H2 

= 

C11H8H+ 

-80.2 

5.6E+14 

34 

H7c9+  +  c2h2 

a 

CHH8H+ 

0.9 

5.6E+14 

35 

H7C9+  +  C2H2 

= 

H/ll+ 

94.1 

5.6E+14 

36 

C9H7+  +  C3H4 

a 

C10H9+  + 

C2H2 

-185.3 

5.9E+14 

37 

c9h7+  +  c3h4 

- 

¥l0+  + 

C2H2 

-54.9 

5.9E+14 

38 

Cgll7+  +  C4H2 

a 

C13H6H3t 

-209.9 

5.8E+14 

39 

H7C9+  +  C4H2 

= 

C13H6H3 

-128.8 

5.8E+14 

40 

C10H9  +  C2H2 

= 

+ 

C12H9  + 

H2 

-23.0 

5.5E+14 

41 

H9C10  +  C2H2 

= 

C12H9  + 

H2 

-153.4 

5.5E+14 

42 

¥l0  +  C2H2 

= 

¥l2+  + 

„2 

-10.5 

5.5E+14 

43 

C10H9  +  C2H2 

= 

H9C12+  + 

Hj 

119.8 

2.7E+06 

7.8E+09 

1.1E+11 

4.2E+11 

44 

H9C11++  M 

= 

Cll¥+  + 

H 

-93.2 

4.5E+14 

45 

C11H8H+  +  C2H2 

a 

C13H6H3 

+  H2 

-122.8 

5.5E+14 

46 

C.1HgH++  C2H2 

= 

C13H8H+ 

+  h2 

86.0 

2.3E+13 

4.0E+12 

2.6E+12 

2.3E+12 

47 

H9C11++  C2H2 

= 

C13H8H+ 

+  H2 

-7.2 

5.5E+14 

5.5E+14 

5.5E+14 

2.8E+14 

48 

H9C11  +  C2H2 

= 

C13H6H3+ 

+  h2 

-216.0 

5.5E+14 

49 

C12H9  +  C2ll2 

= 

C14H8H3+ 

-6.1 

5.5E+14 

50 

¥i2+  +  c2h2 

= 

C14H8H3+ 

-149.0 

5.5E+14 

51 

C12H9++  C4H2 

= 

C16H11+ 

-255.1 

5.6E+14 

52 

H9C12  +  C4H2 

a 

C16H11 

-398.0 

5.6E+14 

53 

¥l2+  +  C4H2 

= 

C16H10H+ 

-135.8 

5.6E+14 

54 

C13H8H+  +  H 

= 

C13H6H3 

+  H 

-208.9 

4.5E+14 

55 

C13H6H3  +  C2H2 

a 

C15H10H+ 

14.2 

5.4E+14 

56 

C13H8H+  +C2H2 

= 

C15H10H+ 

-194.7 

5.4E+14 

57 

C13H8H+  +C4H2 

a 

H11C17+ 

-258.2 

5.6E+14 

58 

C13H8H+  +C4H2 

- 

C17H11+ 

-130.7 

5.6E+14 

59 

C13H6H3+  +  C4H2 

a 

C17H11+ 

78.1 

5.6E+14 

60 

C13H6H3  +  C4H2 

a 

H11C17 

-49.3 

5.6E+14 

61 

C14H8H3  +  C2H2 

= 

C18H11+ 

-242.1 

5.4E+14 
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rcncird  Sate  Cccfficieat.  k 
_ ea^/coTe/s _ 


Mo. 

9  5  A  C  T 

i 

0  S 

A3M7S01 

500 

1CC0 

1500 

2CCC- 

62 

C14H3H3+  4  ^ 

= 

ciSHj(^  +  **2 

20.0 

4.05+14 

2.35*14 

2.25+14 

2.35+14 

63 

C15H10H"  +  C2H2 

= 

C17H11+  4 

«2 

7C.S 

1.45+14 

1.55*13 

3.05*12 

6.15*12 

64 

C15H10H+  +  C2H2 

= 

H11C17+ 

«2 

-55.6 

5.45+14 

65 

C16H11  +  C2H2 

= 

L18H10H*  4  «2 

124.8 

5.45+14 

5.65+13 

6.25*11 

6.35*10 

65 

C16Hir  +C2H2 

= 

H11C18+  +  «2 

112.2 

2.75+11 

2.75*11 

3.45*11 

4.75*11 

67 

C16H10K"+  C2H2 

= 

H11C18+  +  «2 

-149.9 

5.45+14 

S? 

C16H10H++  C2H2 

= 

=  C13H,0H+  +  Hg 

-137.4 

5.45+14 

69 

C16H11++  C4H2 

= 

H11C20++  H2 

-24.4 

5.45+14 

5.45*14 

5.45*14 

5.75*14 

70 

C17H11++  C2H2 

= 

C19H11  4  ^2 

-234.7 

5.35+14 

71 

H11C17  +  C2H2 

= 

C19H11++  H2 

-107.3 

5.35+14 

72 

H11C18++  C2H2 

= 

H11C20^  +  H2 

-129.8 

5.35+14 

73 

C18H10H++  C2H2 

= 

H11C20  4  H2 

-142.4 

5.35+14 

74 

H11C18+  +  C2H2 

= 

C20Hir  +  H2 

50.9 

5.35+14 

5.35+14 

1.05+14 

9.35+12 

75 

C18«10H++  C2H2 

= 

C20H11  +  H2 

36.3 

2.75+13 

4.25+13 

5.85*13 

7.05+13 

76 

Vu*  *  % 

= 

Wl**  "2 

-49.9 

5.35+14 

77 

Vu*  C2B2 

= 

Ws* 

-333.6 

5.35+14 

78 

H11C20  +  C2H2 

= 

C22H13 

-152.9 

5.35+14 

79 

C21H11  +  C2H2 

= 

C23H13 

18.3 

5.35+14 

80 

C21H11+  +  C3H4 

= 

C22H13++  C2H2 

-224.4 

5.55+14 

81 

C22H13  +  C2H2 

= 

C24H13++  H2 

-77.7 

5.35+14 

82 

C22H13  +  C4H2 

= 

C24H13+  +  C2H2 

-84.5 

5.35+14 

c. 

ION-ELECTRON  RECOMBINATION  REACTIONS3 

No. 

REACT 

I  0  11 

A 

n 

1 

H30+  +  e 

- 

H20  +  H 

1.35+19 

-0.5 

2 

HC0+  +  e 

- 

CO  +  H 

7.4E+13 

-0.68 

3 

CH3+  +  e 

- 

CH  +  H2 

5.3E+18 

-0.5 

4 

C2H3+  +  e 

- 

c2h  +  h2 

8.5E+18 

-0.5 

5 

c3H3+  +  e 

-4 

C^  +  CH 

1.1E+19 

-0.5 

6 

H3C3  +  e 

-4 

C2H2  +  CH 

1.1E+19 

-0.5 

7 

C4H5  +  e 

- 

C,H?  +  C,H, 

1.2E+19 

1 

0 

U1 

8 

C^H4  +  e 

•4 

C,H  +  C,H, 

1.3E+19 

-0.5 

9 

h5C5+  +  e 

•4 

C^H^  +  C7H7 

1.4E+19 

-0.5 

10 

C6H5+  *  e 

-4 

c4H4  +  C2H 

1.4E+19 

-0.5 

11 

C?H4H+  +  e 

-4 

C.^Hp  +  C,H, 

1.5E+19 

-0.5 

12 

H7C7+  +  e 

-4 

C6H4  +  CH3 

1.5E+19 

-0.5 

13 

C8H7+  +  e 

-4 

CfiHR  +  C?H 

1.6E+19 

-0.5 

14 

C9H7+  +  e 

-4 

CgHg  +  CH 

1.7E+19 

-0.5 

15 

H7C9+  *  e 

- 

CgHg  +  CH 

1.7E+19 

-0.5 

16 

C10H9+  +  e 

-4 

C10H8  +  H 

1.8E+19 

-0.5 
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ats. 

S  £  A  C 

7  !  0  S 

A 

a 

17 

%CI0+  +  e 

- 

C3^  +  C2H3 

1.22*19 

-0.5 

IS 

CjSH<T*+e 

- 

^"0%  + 

1.52+19 

-0.5 

IS 

*%Fn**e 

- 

1.52*19 

-0.5 

23 

CjjTV  +  e 

- 

Ci2*3  +  H 

1. 52*19 

-0.5 

21 

4C!2*  *  e 

- 

Clfo  *  CjK 

1.52+19 

-0.5 

22 

C13H5r3  *  e 

- 

C13^ 

2.G2+19 

-0.5 

23 

CjjHUHT  +  e 

- 

hrz  *  a 

2.52*19 

-9.5 

24 

C«49ffV  +  e 

- 

C;4H10  +  H 

2.12+19 

-0.5 

25 

cisrHi/  +  e 

- 

C14»10"® 

2.12+19 

-0.5 

25 

+  e 

- 

C16Ki0  +  H 

2.22+19 

-0.5 

27 

cisTi/  +  e 

- 

ClolO  +  H 

2.22+19 

-0.5 

23 

C17Hll^  *  e 

- 

c*5»10  +  ra 

2.22+19 

-0.5 

23 

KilC17+  +  e 

- 

C17H11 

2.22+19 

-0.5 

30 

H.X.*  +  e 

::  io 

- 

C1&H10  +  C2H 

2.32+19 

-0.5 

31 

C13K!0H  *  e 

- 

C13H10  +  K 

2.32+19 

-0.5 

32 

C19H11+  *  e 

- 

C1SH10  +  01 

2.32+19 

-0.5 

33 

c2oHn+  +  e 

- 

C18H10  +  C2H 

2.42+19 

-0.5 

34 

H  C20"  +  e 

- 

C20H11 

2.42+19 

-0.5 

35 

C21Hn+  +  e 

- 

CJ3H10  +  C3H 

2.42+19 

-0.5 

35 

C2ZH13+  +  e 

- 

C22H12  +  H 

2.5E+19 

-0.5 

37 

Vo**' 

- 

c22h,2  +  ch 

2.52+19 

-0.5 

33 

C24H13+  *  e 

- 

Wb*« 

2.52+19 

-0.5 

a  C3K3+  =  linear  H3C3+  =  cyclic  =  allene;  H,Cg  =  propyne;  CgH~+  =  linear  CgH.+; 
linear  C2H3+;  =  benzyl;  =  H-C=C-CH;  CH  =  doublet  delta  state  of  CH  radical;  CxHy+,  KyCx+, 
CxH(y-n)Hn+'  etc‘  rePresent  different  isocers  of  the  sane  ion. 
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nonpolar  species  such  as  the  neutral  growth  species  in  our  ionic  model,  except  for  propyne,  Eq.  (1) 
reduces  to  the  Langevin  equation  which  does  not  have  a  temperature  coefficient.  We  have  used  Eq. 
(1)  to  calculate,  where  necessary,  ion-molecule  reaction  rate  coefficients  for  our  soot  growth  model. 

Several  theoretical  analyses  have  had  the  objective  of  defining  the  temperature  effect  on  ion- 
molecule  reactions  but  these  have  all  concentrated  on  the  situation  when  the  neutral  reactant  has  a 
dipole  moment.17,13,21-24  In  the  Langevin  theory,23  which  forms  the  basis  for  these  analyses,  the  ion 
is  treated  as  a  point  charge.  This  is  certainly  not  'rue  for  many  of  the  large  ions;  even  when  the 
charge  is  localized  it  would  be  shielded  from  the  approaching  reactant  by  the  rest  of  the  molecule. 
Intuitively  the  larger  ions  should  have  a  smaller  rate  coefficient  and  a  negative  temperature 
coefficient;  we  plan  to  improve  the  theoretical  basis  in  an  extension  of  this  program. 

Equation  (1)  accounts  only  for  the  number  of  collisions  (it  does  not  include  collision 
efficiency)  and  fits  most  room  temperature  experimental  rate  coefficient  data.  There  is  considerable 
evidence  that  the  rate  of  ion-molecule  reactions  is  directly  dependent  upon  the  exothermicity  of  the 
reaction,  AHr  and  we  hope  to  use  such  correlations  to  estimate  the  collision  efficiency  for  this  set  of 
reactions  but  have  not  yet  devised  a  logical  means. 

One  of  the  major  problems  in  working  with  large  ions  is  their  identification;  mass 
spectrometry  gives  mass  only.  The  number  of  carbon  and  hydrogen  atoms  has  been  determined  by  the 
use  of  isotopes.3  Thus  for  a  given  molecular  formula  there  can  be  several  isomeric  structures.  We 
thus  include  several  isomers  for  some,  ions  when  their  free  energies  of  formation  are  close;  we  seek 
a  rational  means  of  reducing  this  to  "one  isomer"  per  ion,  probably  by  weighting  the  thermodynamic 
quantities,  AHf  and  Cp  appropriately  and  by  mechanistic  considerations. 

The  present  set  of  ion-molecule  reactions,  chosen  from  a  set  of  about  318  reactions,  is 
presented  in  Table  I.B  and  Fig.  2.  We  have  chosen  to  write  the  reactions  always  toward  increasing 
molecular  size;  thus  the  free  energy  of  reaction  becomes  more  positive  as  the  temperature  is  increased. 
To  avoid  excessively  fast  reverse  reactions,  we  limit  the  forward  reaction  rate  coefficient  so  that  the 
reverse  reaction  never  exceeds  the  Langevin  rate  calculated  by  Eq.  (1).  Because  the  temperature 
coefficient  cannot  be  represented  by  an  Arrhenius  expression,  this  has  caused  some  difficulties  with 
Frenklach’s  code.  Brown  incorporates  the  effect  with  a  table  of  free  energies  as  a  function  of 
temperature. 


C.  ION-ELECTRON  RECOMBINATION 

In  choosing  product  channels  for  the  large  ion-electron  dissociative  recombination  reactions, 
only  molecules  observed  by  Bockhorn  et  al.26  have  been  considered  as  products.  Reaction  rate 
coefficients  for  ion-electron  reactions  are  not  strongly  temperature  dependent,  but  they  do  increase 
with  the  size  of  the  ion. 

The  ions  disappear  by  either  ion-electron  or  positive  ion-negative  ion  recombination  reactions. 
Negative  ion  recombination  rate  coefficients  are  about  one  to  two  orders  of  magnitude  slower  than 
electron  recombination  rate  coefficients.  Further,  negative  ion  concentrations  which  have  been 
measured  are  about  two  orders  of  magnitude  smaller  than  electron  concentrations27  but  there  are  no 
good  measurements  of  electron  or  negative  ion  concentrations  in  soot  forming  flames.  There  is, 
however,  evidence  for  the  presence  of  large  negative  ions  in  soot  forming  flames.5,28,29  We  neglect 
negative  ions  for  the  present;  they  should  be  included  in  the  future.  Their  inclusion  would  decrease 
the  number  of  by-product  neutral  species  which  Frenklach  improperly  uses  to  compare  the  free 
radical  and  ionic  mechanisms.  In  one  computer  run  we  reduced  the  ion  recombination  rate  coefficient 
for  all  reactions  by  two  orders  of  magnitude  and  it  made  a  significant  change  in  the  ion 
concentrations. 
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FIGURE  2  ION  FLOW  DIAGRAM 
Unless  indicated  the  reactant  is  C2H2. 
The  numbers  are  the  free  energy  of  react 
at  1750  K. 
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We  estimate  the  rate  of  ion  recombination,  a,  by  the  equation  for  the  rate  of  collision  of 
electrons  with  particles30 


«d2 

r  ^ 

8kT 

1/2 

r  > 

O 

1  -«■  e" 

4 

'me 

h  j 

(2ire0d)kT 
-  - 

in  which  d  =  the  ion  diameter,  me  =  the  electron  mass,  and  e0  =  dielectric  constant  of  free  space.  The 
ion  diameters  were  calculated  from  ion  mobilities  by  use  of  the  Langevin  equation  for  ion  mobilities. 
Equation  (2)  gives  a  T-1^2  temperature  dependence  which  compares  favorably  with  experiments  of 
Ogram  et  al.31  for  H30+.  The  values  calculated  by  Eq.  (2)  were  about  twice  the  measured  values  of 
Ogram,  so  we  have  divided  the  calculated  values  by  2. 

If  there  is  an  error  in  our  recombination  rates  it  will  be  on  the  high  side.  This  is  consistent 
with  a  comparison  of  experiment  and  computer  runs.  Generally,  in  the  computer  runs  to  date,  most 
large  ions  (there  are  exceptions)  are  computed  to  have  far  too  small  a  concentration  and  a  more  rapid 
decay  rate  than  observed  perimentally,  see  below.  The  small  calculated  concentration  may  be  due 
to  the  ion-molecule  rates  being  too  small,  or  due  to  too  fast  a  loss  rate  by  ion-electron  recombination. 
This  could  indicate  the  presence  of  significant  concentrations  of  negative  ions;  their  rates  of 
recombination  are  smaller  than  for  electrons.  Of  course,  the  recombination  of  a  large  positive  ion  and 
a  large  negative  ion  would  double  the  size  of  the  carbon  species,  albeit  neutral,  and  could  be 
interpreted  as  the  first  step  in  coagulation.  Such  collisions  are  favored,  by  electric  charge  effects, 
over  ion-neutral  or  neutral-neutral  collisions. 


V.  DIFFUSION 


Because  the  pressure  for  the  standard  flame  is  less  than  one  atmosphere  and  we  have 
experimentally  demonstrated  the  importance  of  diffusion5  in  the  flame  chosen  for  study,  it  is 
necessary  to  include  the  diffusion  coefficients  for  ions.  We  estimated  these  using  the  procedure 
previously  developed  for  interpreting  Langmuir  probe  data.7  Experimental  ion  mobilities  of  a  wide 
mass  range  of  PCAH  ions32,33  were  extrapolated  to  the  higher  mass  range  required  for  this  program. 
This  extrapolation  gives  results  of  /r  vs.  ionic  mass,  amu,  which  parallels  the  results  obtained  from  a 
calculation  using  the  Langevin  type  equation,  giving  confidence  in  the  procedure. 

We  assume  that  the  negative  species  is  a  free  electron,  although  there  is  some  evidence  that 
there  are  large  negative  ions  present.6,28  With  electrons  present  the  ambipolar  diffusion  coefficient 
of  the  ion  must  be  used  because  the  ions  do  not  diffuse  independent  of  electrons.  The  ambipolar 
diffusion  coefficient  (Da)  for  a  specific  ion  was  calculated  from  the  ion  mobility  using  the 
relationship34:  Da  =  2kT/x/e,  where  k  is  Boltzmann’s  constant,  T  is  temperature,  m  is  the  ion  mobility, 
and  e  is  the  elementary  charge.  The  ionic  mobility  is  calculated  by  the  following  procedure: 

(1)  Calculate  the  ion  mobility,  ja,  at  2.67  kPa  and  273  K  for  each  of  the  dominant  flame  gases 
using  the  following  correlations  developed  at  AeroChem  based  on  the  Langevin  equation  combined 
with  experimental  values: 


HW}  (P0,  T0)  -  a(j)  x  (MW;)b(j)  cm2  V'1  s"1  (3) 
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where  i  is  the  ion  species  of  mass  MWj  and  where  a(j)  and  b(j),  for  each  flame  gas,  j,  are  given  in 
the  following  table:  r2  is  the  coefficient  of  determination  of  the  fit  to  the  actual  data: 


Gas 

a 

b 

r2 

h2 

3471.53 

-0.49893 

0.9981 

o2 

742.83 

-0.48710 

0.9981 

C02 

493.44 

-0.46446 

0.9975 

C2H2 

473.39 

-0.43385 

0.9991 

CO 

688.39 

-0.47242 

0.9989 

h2o 

760.55 

-0.46498 

0.9984 

(2)  Correct  each  of  the  low 

pressure,  low  temperature  ion  mobilities  calculated  above,  to  flame 

pressures  and  temperatures,  by  the  following  formula  (except  for  water),  where  P  is  the  flame 
pressure  (kPa)  and  T  is  temperature  of  the  flame: 


ft  j  (P,  T) 


*  ftj  (po>To) 


(4) 


For  water  vapor: 


(p>  T) 


1.16 


x 


/*i,H20  (P0>  T0) 


(5) 


(3)  Correct  the  ion  mobility  calculated  for  each  of  the  pure  gases,  and  evaluated  at  flame 
conditions,  to  the  actual  flame  gas  composition  using  Blanc’s  Law  -  i.e.,  a  simple  viscosity  mixing 
rule; 


/h"1  (for  the  gas  mixture)  -  Y'  — i  (6) 

i  ftj 


where  Xj  is  the  mole  fraction  of  gas  j  in  the  flame,  and  ju,  is  the  overall  ion  mobility  of  the  ith  ion 
in  the  gas  mixture. 

It  was  apparently  difficult  to  incorporate  this  into  Frenklach’s  computer  program  so  that  he 
has  employed  a  different  procedure  for  including  the  diffusion  coefficient.4 


VI.  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  ION  CONCENTRATION  PROFILES 


Frenklach  and  Wang4  have  run  a  limited  number  of  computer  runs  using  the  sets  of  reactions 
we  developed  for  the  ionic  mechanism  of  soot  formation.  The  early  runs  showed  great  differences 
between  experimental  and  computed  ion  profiles.  As  the  model  was  improved  the  agreement 
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improved.  There  were  no  arbitrary  alterations  of  rate  constants  or  of  thermodynamics.  The  major 
changes  were  in  the  choice  of  reactions  and  in  the  isomer  used;  large  ions  have  a  large  number  of 
isomers,  so  ...ns  gives  a  large  range  of  choices.  Without  making  such  choices  the  number  of  reactions 
and  reactants  becomes  excessive  for  the  computer  memory.  The  effect  of  such  choices  on  the  realism 
of  models  has  yet  to  be  evaluated. 

Some  selected  comparisons  between  calculated  and  experimental  ion  profiles  are  presented  in 
Figs.  3  and  4.  For  small  ions,  Fig.  3,  the  agreement  is  reasonably  good.  The  relationship  of  the  two 
isomers  H3C3+  and  C3H3+,  the  cyclic  and  linear  isomers  of  the  proposed  precursor  ion,  demonstrates 
that  equilibrium  maintains  a  relatively  high  concentration  of  the  reactive  isomer.  One  of  the 
arguments  against  the  ionic  mechanism  was  that  the  stable  isomer  H3C3+  was  found  experimentally 
to  be  nonreactive  with  acetylene  or  other  small  hydrocarbons35,36  and  thus  because  of  this  ion’s 
stability  ionic  growth  would  not  proceed.  The  linear  isomer  was  found  to  be  very  reactive.35,36 

The  data  for  larger  ions.  Fig.  4,  are  not  in  such  good  agreement.  The  peak  concentrations  are 
reasonably  close  but  the  profiles  differ  markedly.  This  difference  is,  however,  no  greater  than  the 
difference  for  neutral  species  in  the  free  radical  mechanism,  where  it  is  declared  as  good  agreement. 
One  major  difference  between  our  approach  to  modeling  and  Frenklach’s  is  that  we  have  limited  the 
reactants  and  products  to  those  that  have  been  observed  experimentally  and  have  required  that  the 
model  calculate  all  such  species.  Frenklach’s  model  does  not  have  either  of  these  constraints.4  Thus 
his  model  includes  cyclopentaphenanthrene  which  is  not  observed  experimentally,  and  does  not 
include  phenanthrene  and  pyrene  which  are  observed  experimentally.  Further,  in  Frenklach’s  model,4 
only  five  calculated  cyclic  compound  profiles  are  compared  with  experiment  and  the  peak 
concentration  of  one  of  these,  naphthalene,  is  computed  to  be  20  times  smaller  than  measured. 

The  calculated  ion  concentrations  increase  with  distance  and  the  experimental  concentrations 
decrease,  Fig.  4.  A  calculation  was  carried  out  in  which  the  temperature  was  fixed  at  the  maximum 
temperature  (the  temperature  profile  is  an  input  parameter);  the  calculated  concentrations  then 
decayed  slowly.  This  shows  the  sensitivity  to  temperature.  In  an  early  computer  experiment  a  point 
was  misplaced  in  the  temperature  rise  section  of  the  temperature  profile  and  this  produced  a 
corresponding  change  in  the  ion  profile,  confirming  the  temperature  sensitivity  of  the  computer 
model. 


We  compare,  in  Fig.  5,  the  maximum  concentration  of  ions  experimentally  observed  and 
calculated.  Since  the  experiments  do  not  distinguish  between  isomers  and  sometimes  there  are  more 
than  one  isomer  in  the  computer  program,  we  compare  in  Fig.  5  the  isomer  in  which  the  agreement 
is  best  between  the  calculated  and  experimental  concentrations.  For  most  of  the  larger  ions,  the 
experimental  value  exceeds  the  calculated  value.  We  have  neglected  electron  attachment  and  thus 
positive  ion-negative  ion  recombination,  which  if  included  would  increase  the  calculated 
concentrations;  ions  would  not  be  removed  from  the  system  as  fast  as  by  ion-electron  recombination. 
Incidentally  this  would  decrease  the  formation  of  neutral  species  which  Frenklach  incorrectly 
interprets  as  a  measure  of  the  efficacy  of  the  ion-molecule  mechanism. 

The  free  radical  mechanism  which  is  used  to  produce  the  precursor  ion  in  the  computer  model 
does  not  agree  well  with  experimental  measurements  of  these  neutral  species,  Fig.  6.  The 
experimental  measurements  have  been  made  by  several  people  with  good  agreement.  The  question 
is  thus  raised;  does  the  difference  between  experiment  and  calculation  for  the  ionic  mechanism 
depend  upon  the  problems  in  the  free  radical  mechanism  or  are  they  problems  in  the  ionic 
mechanism?  This  could  be  readily  tested  if  the  experimental  neutral  species  profiles  required  for  the 
ionic  mechanism  could  be  an  input  to  the  computer  program,  just  as  the  temperature  is.  Frenklach 
could  not  handle  this  with  his  program  so  the  test  was  done  by  Brown  and  Pedersen  whose  computer 
program  does  have  this  capability.  First  they  ran  the  program  using  as  input  the  calculated  neutral 
species  concentrations  from  Frenklach.  The  resulting  ion  profiles  were  very  close  to  those  calculated 
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by  Frenklach,  confirming  the  validity  of  the  technique  of  using  concentration  profiles  as  input  to  the 
computer  program. 

When  Brown  and  Pedersen  then  ran  the  program  using  as  input  the  experimental  neutral 
species  profiles  required  by  the  ionic  mechanism,  the  maximum  calculated  ion  concentrations  were 
many  orders  of  magnitude  lower  than  measured!  This  indicates  a  great  sensitivity  in  the  computer 
model  to  the  free  radical  mechanism.  One  has  to  ask  the  question:  how  meaningful  are  such  models 
if  they  are  so  sensitive?  This  question  has  yet  to  be  answered;  more  sensitivity  analyses  are  clearly 
called  for  and  will  be  done  in  the  extension  of  the  present  work.  The  only  real  way  to  test  a  model 
is  to  apply  it  to  more  than  one  system.  We  also  plan  to  do  this  in  the  extension  of  this  work  when  we 
will  apply  this  model  to  the  benzene  flame. 


VII.  COMPARISON  OF  THE  RATE  OF  CARBON  SPECIES  GROWTH 
FOR  THE  FREE  RADICAL  AND  IONIC  MECHANISMS 


A  simple  technique  has  been  employed  to  compare  the  relative  rates  by  which  the  two 
competing  mechanisms  account  for  the  growth  of  carbon  containing  species  in  the  same  sooting  flame. 
It  is  accepted  that  soot  is  formed  from  small  carbonaceous  molecules  by  a  mechanism  in  which  they 
increase  in  size  and  carbon  to  hydrogen  ratio  until  at  some  size  they  coagulate  to  produce  even  larger 
species,  and  eventually  produce  incipient  soot  particles.  The  least  understood  steps  in  this  process  are 
the  increase  in  carbon  number  from  two  to  a  molecular  weight  of  about  1,000  u.  In  this  comparison 
we  use  experimental  data  combined  with  reaction  rate  coefficients  to  compare  the  rates,  and  thus 
times,  for  the  addition  of  ten  carbon  atoms  to  the  growing  molecular  species  in  the  range  of  carbon 
numbers  for  which  experimental  data  are  available.  The  range  of  observed  carbon  numbers  is  very 
limited  for  neutral  species  because  as  they  become  larger  their  concentration  rapidly  diminishes  below 
detectability. 

For  this  comparison  we  chose  the  well  studied  acetylene/oxygen  flat  flame  at  2.67  kPa, 
equivalence  ratio  =  3.0  and  unburned  gas  velocity  =  50  cm/s.  For  the  experimental  neutral  species 
concentrations  we  use  those  of  Vovelle3  because  they  extend  to  mass  252  u.  Bittner  and  Howard’s 
data8  extend  only  to  mass  178,  C14H10,  and  at  mass  178  are  more  than  an  order  of  magnitude  less  than 
those  of  Vovelle.  This  choice  favors  the  free  radical  mechanism.  For  the  free  radical  mechanism, 
free  radical  concentrations  have  not  been  measured,  so  they  are  assumed  equal  to  the  concentration 
of  the  preceding  neutral  species  from  which  the  radical  was  formed.  This  again  favors  the  free 
radical  mechanism. 

For  the  individual  ionic  species  concentrations,  we  use  those  measured  at  AeroChem5;  the  total 
concentrations  have  been  confirmed  by  Gerhardt  and  Homann29  and  are  consistent  with  measurements 
of  Delfau,  Michaud  and  Barassin.37 

For  the  free  radical  mechanism  we  use  that  of  Frenklach  and  Warnatz38  and  use  the  rate 
coefficients  of  Frenklach  and  Wang39  for  the  specific  steps.  The  free  radical  mechanism  involves  four 
elementary  reactions  with  the  following  rate  coefficients: 


R  is  the  large  reactant  species. 

R  +  H.  (-H,) 

R  +  C2H2  (-  H.) 

R  +  C2H2 
R  +  H. 


k 

k 

k 

k 


4.5 

X 

1012  cm3  mol' 

3.1 

X 

1012 

1.0 

X 

1013 

1.3 

X 

1013 
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The  ionic  mechanism  for  growth  requires  two  reactions: 

R+  +  C2H2(-H2)  k  =  3.6  -  5.3  x  1014 

R+  +  C2H2  k  =  5.6  x  1014 

The  ion-molecule  reaction  rate  coefficients  are  calculated  by  Langevin  theory  for  each 
reaction  and  are  adjusted  so  that  the  reverse  reaction  rate  calculated  by  thermodynamic  equilibrium 
never  exceeds  the  Langevin  rate  calculated  for  the  reverse  reaction. 

The  times  for  the  species  to  add  a  specified  number  of  carbon  atoms  by  the  two  mechanisms 
were  then  compared.  Thus,  for  the  reaction: 

A  +  B  C  +  D 

the  rate  of  reaction  is  given  by: 

R  -  —  m  kAB  (7) 

dt 

For  this  analysis  we  use  the  maximum  (with  respect  to  distance  from  the  burner)  experimental 
concentrations  of  A,  and  use  the  measured  value  of  B  at  the  position  in  the  flame  where  the 
concentration  of  A  is  maximum.  The  appropriate  rate  coefficient  for  k  is  used  for  either  the  free 
radical  or  the  ionic  mechanism.  In  the  free  radical  mechanism,38,40  A  is  a  stable  species  or  a  free 
radical  and  C  is  a  free  radical  with  the  same  number  of  carbon  atoms  as  A,  or  a  stable  species  with 
two  more  carbon  atoms  than  A;  B  is  a  hydrogen  atom  or  acetylene.  In  the  ionic  mechanism,  A  is  an 
ion  and  C  is  an  ion  with  two  more  carbon  atoms  than  A;  B  is  acetylene. 

The  reaction  time  is  calculated  from  the  experimental  maximum  concentration  of  the  growing 
species,  the  concentration  of  the  smaller  reactant,  H*,  H2,  or  C2H2,  at  the  distance  in  the  flame  at 
which  the  growing  species  maximizes  and  the  reaction  rate  coefficient,  by: 

r  =  n/R  (8) 

where  n  is  the  number  of  reacting  species  flowing  through  the  system,  i.e.,  4  x  1014  cm-3;  and  R  is 
the  reaction  rate  for  any  step. 

The  times  required  for  4  x  1014  cm-3  growing  neutral  or  ionic  species  to  add  ten  carbon  atoms 
are  displayed  in  Fig.  7  for  each  step  in  the  mechanism  and  for  the  total  time.  The  time  for  the  reverse 
reaction  for  each  step  is  also  given.  The  reverse  time  is  the  rate  of  the  reverse  reaction  calculated 
from  thermodynamic  equilibrium.  The  number  density  of  reactive  species  of  any  given  variety,  4  x 
1014  cm'3,  is  chosen  because  this  is  the  maximum  number  of  soot  particles  observed  in  this  flame. 

The  total  times  to  add  ten  carbon  atoms  by  the  two  mechanisms  are  comparable:  8.1  /rs  for  the 
free  radical  mechanism  and  6.7  ns  for  the  ionic  mechanism. 

These  total  times  do  not  include  reverse  reactions  which  should  be  taken  into  account;  a  rapid 
reverse  reaction,  small  time,  can  effectively  reduce  the  forward  rate  of  reaction,  depending  on  the 
rates  and  equilibria  of  the  preceding  and  following  reaction  steps.  The  reverse  reactions,  given  in  Fig. 
7,  appear  to  be  a  greater  complication  in  the  free  radical  mechanism  than  in  the  ionic  mechanism. 
An  indication  of  the  forward  reaction  rate  reduction  is  given  by  the  ratio  of  forward  to  reverse 
reaction  times,  when  this  value  exceeds  1,  the  reaction  proceeds  more  rapidly  in  the  reverse  reaction 
than  in  the  forward  direction,  of  course,  the  overall  effect  is  more  important  when  the  forward 
reaction  time  for  any  individual  step  is  long  compared  to  the  average  reaction  time. 
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The  above  analysis  does  not  permit  a  clear  decision  between  the  free  radical  mechanism  and 
the  ionic  mechanism.  For  modest  size  carbon  species,  10  to  20  carbon  atoms,  the  rate  of  growth  of 
carbon  species  is  about  the  same  for  the  free  radical  mechanism  and  for  the  ionic  mechanism.  The 
greater  concentration  of  neutral  species  is  balanced  by  the  greater  reaction  rate  coefficients  for  ion- 
molecule  reactions  and  the  fewer  number  of  steps  involved  in  adding  a  specific  number  of  carbon 
atoms  to  the  growing  species  for  the  ionic  mechanism,  five  steps,  than  for  the  free  radical  mechanism, 
eleven  steps.  The  free  radical  mechanism  suffers  in  the  comparison  because  of  lack  of  experimental 
data  for  the  free  radicals  involved  and  for  larger  molecular  species.  The  free  radical  mechanism  also 
appears  to  have  greater  complications  than  the  ionic  mechanism  due  to  reverse  or  equilibrium 
reactions. 

Consider  the  implications  of  the  above  calculation  of  reaction  times  to  determine  the 
reasonableness  of  this  approach.  The  maximum  soot  number  density,  Fig.  2  in  Ref.  9,  is  reached  at 
about  35  mm  above  the  burner,  or  about  6.7  ms  from  the  position  in  the  flame  at  which  large  carbon 
containing  species  maximize.  This  is  a  good  estimate  of  the  time  available,  rg,  for  soot  particles  to 
be  formed  from  molecular  species.  If  we  assume  the  time  for  the  addition  of  one  carbon  atom  to  the 
growing  species  is  rc,  then  the  number  of  carbon  atoms,  Nc,  that  can  be  added  to  the  growing  nuclei 
is: 


6.7  x  IQ'3 
7.4  x  10-7 


9,000  carbon  atoms. 


(9) 


rc  is  taken  as  the  average  for  the  neutral  and  ion  mechanisms.  9,000  carbon  atoms  corresponds  to  a 
molecular  weight  of  about  110,000  u,  equivalent  to  a  particle  diameter  of  about  4.5  or  3.0  nm, 
depending  upon  whether  the  particle  is  planar  or  spherical,  respectively.41  The  experimentally 
observed  particle  diameters  at  35  mm  above  the  burner  surface  are  9-13  nm  for  neutral  particles  and 
3-6  nm  for  charged  particles.  It  is  interesting  that  the  calculated  diameter,  assuming  the  equivalent 
of  a  fixed  rate  (fixed  time)  for  adding  carbon  atoms  to  the  growing  species,  neutral  or  ion,  leads  to 
a  diameter  of  the  carbon  particle  very  close  to  that  observed,  within  the  accuracy  of  the  calculation 
and  the  measurement. 

The  above  discussion  demonstrates  that  examining  relative  growth  of  small  carbon  species  does 
not  permit  a  conclusion  concerning  the  rate  of  soot  nucleation  by  the  free  radical  and  the  ionic 
mechanism;  they  each  have  about  the  same  reaction  times  using  the  available  data.  To  make  a 
distinction  between  the  two  mechanisms  one  must  examine  reactions  involving  larger  numbers  of 
carbon  atoms. 

Note  the  decrease  in  maximum  concentration  in  going  from  C6H6  to  C14H8  in  the  standard 

flame: 


benzene 

c6h6 

7.8  x  1012 

naphthalene 

c10h8 

1.8  x  1011 

ethynylnaphthalene 

*-"12^8 

6.9  x  1010 

Ci4H8 

1.7  x  1010 

This  is  a  decrease  in  concentration  of  about  a  factor  of  60  per  carbon  atom  added.  Bockhorn’s 
results,26  which  cover  a  larger  range  of  carbon  species  do  not  indicate  such  a  rapid  decay,  but 
unfortunately  they  are  in  a  different  flame. 

For  ions,  the  decay  in  concentration  is  much  less,5  e.g.: 
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C3H3*  4.7  x  10s  per  cm  3 

CjftHts^  1-9  x  107 

c45h;/  1.9  xio6 

This  is  a  decrease  in  concentration  of  about  a  factor  of  6  per  carbon  atom  added,  about  one-tenth  that 

indicated  above  for  free  radical  species! 

Clearly  a  comparison  of  the  two  mechanisms  requires  extension  to  much  larger  species  than 
is  done  here.  The  computer  modeling  effort  will  thus  have  to  be  extended  to  greater  masses  than  arc 
currently  in  the  computer  model. 


VII!.  EXPERIMENTAL  RATES  OF  ION  FORMATION 


One  requirement  of  any  mechanism  of  soot  formation  is  that  the  rate  of  formation  of  soot 
precursors  equal  or  exceed  the  rate  at  which  soot  is  produced.  In  this  section  we  examine  previously 
obtained3  experimental  data  for  the  total  ions  and  the  individual  ions. 

First  we  examine  the  experimental  data  presented  in  Fig.  8  to  compare  the  experimental  rate 
at  which  the  total  number  of  ions  is  generated  with  the  rate  at  which  neutral  soot  is  observed  to  be 
formed.  The  temperature,  total  ion  concentration,  and  the  neutral  and  charged  soot  particle 
concentration  profiles  are  presented  in  Fig.  8.  Identifiable  soot  particles,  (i.e  those  which  can  be 
detected  using  an  electron  microscope,  diameter  exceeding  about  1.5  nm)  first  appear  at  about  2.0  cm 
above  the  burner,  yet  a  yellow  glow,  presumably  due  to  soot,  first  appears  at  about  1.0  cm.  To  reflect 
this  observation,  i.e.,  that  soot  is  formed  where  the  yellow  glow  appears  even  though  the  particles  are 
too  small  to  be  detected  with  the  electron  microscope„we  have  drawn,  in  Fig.  8,  an  interpolated 
(dotted)  soot  concentration  curve  starting  at  1  cm  and  extending  to  the  measured  maximum. 

We  treat  the  flame  as  a  steady  state,  one-dimensional  system.  The  continuity  equation 
describing  the  ion  concentration  at  any  distance  from  the  burner  is: 

(dl/dt)  =  (net  ion  production  rate)  -  V(dl/dx)  +  D(d2I/dx2)  =  0  (10) 

where  I  =  ion  concentration,  V  =  flow  velocity,  D  =  ion  diffusion  coefficient,  and  x  is  the  distance 
from  the  burner.  V  has  been  determined  as  a  function  of  distance  in  this  flame3  and  D  was  calculated 
from  estimated  ion  mobilities,  p,6  using  the  Einstein  relation:  D  =  p(kT/e),  where  k  is  the  Boltzmann 
constant  and  e  the  electronic  charge.  Combining  these  values  with  the  ion  concentration  derivatives 
of  the  profile  in  Fig.  8  gives  the  "net  ion  production  rate"  as  shown  in  Fig.  9. 

The  total  ion  production  rate,  q,  was  obtained  by  adding  the  calculated  ion  loss  rate  by 
recombination  with  f ree  electrons,  to  the  "net  ion  production  rate”.  The  ion  lecombination  rate  is  cd2, 
where  a  is  the  ion-electron  recombination  coefficient.  For  a  we  used  2  x  10"7  cm3  s-1  where  small 
(i.e.,  39  amu)  ions  dominate,  and  corrected  for  increasing  ion  mass  downstream  in  the  flame6  using 
a  factor  proportional  to  d2  where  d  =  ion  diameter  and  a  capacitive  term  of  the  form  (1  +  A/d)  with 
A  a  constant.  The  total  ion  production  rate,  q,  is  also  plotted  in  Fig.  9  as  a  dashed  line— note  scale 
change.  If  negative  ions  are  assumed  to  be  the  recombining  partners  for  positive  ions,  the  total  ion 
production  rates  would  be  reduced  because  ion-ion  recombination  is  slower  by  almost  an  order  of 
magnitude  than  ion-electron  recombination. 

The  net  ion  production  rate  shows  two  peaks  corresponding  to  the  two  peaks  in  the  ion 
concentration  curve,  the  source  of  these  two  peaks  is  unknown.  Also  shown  in  Fig.  9  is  the  net  rate 
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of  particle  formation  derived  from  the  two  neutral  soot  curves  in  Fig.  8  using  Eq.  (10),  in  which  the 
diffusion  term  for  these  large  particles  is  now  negligible  compared  with  the  flow  velocity  term. 

The  first  observation  is  that  the  maximum  net  rate  of  ion  f ormation,  1 .5  x  1 013  ions  cm-3  s— \ 
exceeds  the  net  rates  of  soot  formation,  2.5  and  7.5  x  1012  particles  cm-3  s— 1  from  either  curve. 
Second,  the  net  ion  loss  and  soot  formation  occur  in  the  same  region  of  the  flame.  Third,  the  net  ion 
loss  rate,  about  2  x  1012  cm-3  s- \  corresponds  with  the  net  particle  formation  rate,  about  2  x  1012 
cm-3  s-1  assuming  particles  first  appear  at  the  position  in  the  flame  where  yellow  first  occurs  (dotted 
curve.  Figs.  8  and  9).  The  peak  particle  production  rate  for  the  measured  soot  curve  is  only  about 
four  times  greater  than  the  peak  ion  net  loss  rate.  The  value  for  the  net  ion  disappearance  rate, 
however,  does  not  include  larger  ions  lost  by  recombination  which  can  still  lead  to  soot  formation  via 
the  same  types  of  reactions  as  in  the  postulated  free  radical  mechanisms,  see  Fig.  1 . 

These  analyses  lend  further  support  to  an  ionic  mechanism  of  soot  formation  in  this  flame. 
Clearly  more  precise  data  are  desirable,  particularly  on  negative  molecular  ions  and  on  their  rates  of 
recombination,  and  particle  concentration  and  particle  size  distribution  (both  neutral  and  charged) 
where  soot  is  first  observed,  i.e.,  where  the  yellow  emission  first  appears. 

The  experimental  rates  of  formation  and  destruction  for  three  ions,  C3H3+,  C13H9+,  and 
C21Hn+  are  shown  in  Figs.  10-12.  In  all  three  examples,  the  maximum  rates  of  formation  and 
destruction  are  the  order  of  1012  ions  cm-3  s— \  consistent,  within  the  experimental  accuracy,  with 
the  rate  of  soot  formation.  Further,  these  small  ions  start  to  disappear  almost  exactly  at  1.0  cm  where 
soot  is  first  observed  in  this  flame.  These  observations  are  again  consistent  with  the  ionic  mechanism 
of  soot  formation.  A  similar  comparison  cannot  be  made  for  the  free  radical  mechanism  because 
measurements  are  not  available. 


IX.  RECOMMENDATIONS  FOR  FUTURE  WORK 


Future  work  on  soot  formation  should  be  concentrated  on  solving  problems  as  opposed  to 
producing  more  measurements.  Computer  modeling  is  the  basis  for  interpreting  the  mechanism,  from 
nucleation  to  final  agglomeration;  experimental  work  should  be  driven  by  the  need  of  modelers. 
Understanding  the  nucleation  step  remains  the  major  challenge,  all  of  the  other  steps  in  the  process 
are  reasonably  well  understood. 

The  major  experimental  need  for  the  free  radical  mechanism  is  for  information  on  profiles 
of  species  larger  than  mass  252  u,  five  aromatic  rings.  This  is  the  largest  neutral  species  profile 
available  in  a  sooting  flame  to  compare  with  models.  If,  as  appears  to  be  the  case,  the  experimental 
measurements  are  unavailable  because  the  concentrations  become  too  small  to  measure  when  the  size 
exceeds  252,  then  there  is  reason  to  question  the  validity  of  the  free  radical  mechanism.  Reaction 
rates,  for  such  small  concentrations,  would  be  too  small  to  produce  soot  at  the  observed  rates  of 
production.  This  is  not  an  easy  measurement  to  make.  It  will  have  to  be  done  with  a  molecular  beam 
sampling  mass  spectrometer  which  uses  photoionization.  Commonly  used  thermal  electron 
fragmentation  of  the  parent  species  makes  it  difficult  to  identify  the  parent  species.  Both  stable  and 
free  radical  species  must  be  measured. 

Mass  spectrometry  is  an  intrusive  technique  and  has  thus  been  n^6lected,  the  major  effort  has 
been  on  nonintrusive  techniques.  Expensive  and  time  intensive  laser  techniques  have  not  been 
capable  of  identifying  and  quantifying  the  PCAII  species  in  sooting  flames.  The  spectra  obtained  are 
too  diffuse  to  allow  this  type  of  analysis.  It  is  now  time  to  use  the  best  available  experimental 
techniques  to  obtain  the  required  data. 
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FIGURE  12  NET  EXPERIMENTAL  PRODUCTION  RATE  OF  C21Hn+ 


Another  major  experimental  need  which  is  not  being  addressed  is  the  rate  at  which  incipient 
soot  particles  are  formed.  The  concentration  profiles  of  species  falling  between  large  ions  (and  large 
neutral  species)  and  small  soot  particles  are  needed.  This  is  also  difficult  to  measure.  The  data  are 
however  very  important  for  understanding  both  the  free  radical  and  the  ionic  mechanisms. 

The  major  experimental  need  for  the  ionic  mechanism  is  concentration  profiles  of  negative 
ions  in  both  the  standard  acetylene  flame  and  in  the  benzene  flame  in  which  Howard  and  associates 
have  measured  neutral  species  profiles  and  AeroChem  has  measured  positive  ion  species  profiles. 

For  all  experimental  work,  it  is  extremely  important  that  a  variety  of  measurements  be  made 
on  the  same  flame  by  a  number  of  researchers  to  reduce  the  degree  of  freedom  which  modelers 
currently  enjoy  because  of  limited  data  on  any  given  system.  Unless  a  flame  is  well  documented, 
which  no  one  research  group  can  do,  experimental  measurements  are  essentially  worthless.  A  few 
experimental  measurements  on  a  unique  flame  can  be  used  to  prove  most  anything  and  can  very  easily 
be  simulated  by  models.  The  conclusions  are  worthless.  Unfortunately  the  literature  is  replete  with 
such  measurements  and  interpretations. 
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APPENDIX  A 


A  Computer  Code  Designed  for  Modeling  the  Ionic 
Mechanism  of  Soot  Formation 

Robert  C.  Brown  and  Timothy  W.  Pedersen 
Iowa  State  University 

As  a  prerequisite  for  this  study,  it  was  necessary  to  develop  a  versatile  computer 
code  that  could  incorporate  features  that  are  not  standard  in  reactive  flow  models. 
In  particular,  the  code  developed  for  this  study  can  treat  non-Arrhenius  kinetics 
by  use  of  a  lookup  table.  The  code  also  treats  ionic  mobility  by  calculating  the 
electric  field  distribution  throughout  the  flame  region.  Changes  in  reaction  rates  or 
reaction  mechanisms  were  easily  implemented  into  the  computer  code  because  of  this 
versatility. 

The  reaction  mechanism  used  in  this  study  was  separated  into  two  models:  a 
neutral  species  model  and  an  ionic  species  model.  This  partition  is  based  on  the 
assumption  that  the  ionic  reactions  do  not  change  the  neutral  species  and  temperature 
profiles  significantly.  The  neutral  species  profiles  are  treated  as  parameters  in  the 
ionic  model. 

Comparison  of  predicted  and  experimental  neutral  species  profiles  for  rich  acety¬ 
lene  flames  demonstrated  significant  deficiencies  in  the  present  understanding  of  neu- 


32 


tral  species  mechanisms.  Since  this  study  focuses  on  ionic  mechanisms  of  soot  nucle- 
ation,  we  chose  to  use  experimental  profiles  of  neutral  species  as  input  to  the  ionic 
model.  Furthermore,  since  the  cyclic  form  of  C3H3  is  the  starting  ion  for  the  pro¬ 
posed  ionic  mechanism,  the  experimental  profile  of  this  ion  was  also  used  as  an  input 
parameter  to  the  ionic  model,  thereby  circumventing  the  need  for  chemiionization 
reactions  in  the  ionic  model. 

I.  GOVERNING  EQUATIONS 

The  equations  describing  a  one-dimensional  laminar  premixed  flame  in  the  pres¬ 
ence  of  an  electric  field  include  the  mass  and  species  conservation  equations  and 
Poisson’s  equation  to  describe  the  electric  field  intensity. 


Mass  conservation: 


dp  ■  d(pu) 

dt  dx 


=  0 


1  d{AJi )  , 
~A~dT  +  U' 


(A.1) 

(A.2) 


where  Y;  is  the  specie  mass  fraction  and  ux  is  the  rate  of  chemical  generation.  The 
mass  flux  J;  is  defined  as 


n  y  ■ 

Ji  =  +  KmpYiE.  (A.3) 

The  first  term  on  the  right  side  of  the  equation  is  the  mass  flux  due  to  molecular 
diffusion  and  the  second  term  is  the  mass  flux  due  to  ion  mobility  in  the  electric  field. 
To  avoid  the  use  of  ambipolar  diffusion  coefficients,  a  new  method  for  calculating  the 
mobility  of  ions  in  the  presence  of  a  self-induced  electric  field  is  used  in  this  study.  To 


33 


calculate  the  mass  flux  due  to  ion  mobility,  the  electric  field,  E,  is  calculated  directly 


from  Poisson’s  equation: 


<PV  1 

—  =  -_(n+-n.)e 


(A-4) 


where  V  is  the  voltage,  e0  is  the  permissivity  constant,  e  is  the  electron  charge,  and 
n+  and  n_  are  positive  and  negative  ion  concentrations,  respectively. 

The  gas  velocities  for  the  flames  in  this  research  are  on  the  order  of  100  cm/s. 
The  pressure  throughout  the  flame  is  assumed  to  be  nearly  constant.  Because  of  the 
small  velocity  and  small  pressure  drop,  the  momentum  equation  does  not  need  to  be 
solved.  Since  an  experimental  temperature  profile  was  used,  there  was  no  need  to 
solve  the  energy  equation. 

To  eliminate  the  convection  terms  from  the  conservation  equations  and  to  auto¬ 
matically  satisfy  continuity,  a  coordinate  transformation  [2]  is  made  using 


dip  ,  dip 

ani  a =  ~pu' 


The  transformed  equation  for  the  ionic  species  mechanism  is 


(A.5) 


dy 

'it 


A  (nV^t 

dip  V  dip, 


P~KT  ~  PlTJ.  \  P^-ZJ.  +  P7u.W)  +  G 


d 

dip' 


(A.6) 


where  the  expressions  corresponding  to  T,  T',  and  G  are  in  Table  A.l.  The  appropriate 
boundary  conditions  are 

</3  =  V?o  at  ip  —  0  (A.7) 


and 


d<p 

dip 


—  constant  at  ip  =  oo. 


(A.8) 


The  variable  <p  is  a  generic  term  for  Y;  and  ip  is  the  transformed  spatial  variable. 
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Table  A.l:  Values  of  V,  V  and  G 
for  the  ionic  species 
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Equation  (A.  10)  represents  the  time  rate  of  change  of  <p  due  to  transport  and  Equation 
(A.li)  is  the  time  rate  of  change  of  e?  due  to  chemical  generation. 

Equations  (A.10)  and  (A.li)  are  partial  differential  equations  (PDEs)  in  time 
and  space.  Partial  differential  equations  are  by  nature  more  difncult  to  solve  than 
ordinary  differential  equations  (ODEs).  A  system  of  ODEs  would  have  the  added 
advantage  of  using  an  ordinary  differential  equation  solver  in  the  solution  process. 
The  method  of  lines  converts  the  system  of  PDEs  to  a  system  of  ODEs  [4].  This 
technique  discretizes  the  spatial  derivatives  leaving  a  system  of  first  order  ODEs  in 
time  [2].  Central  differencing  of  Equations  (A.lO)'and  (A.li)  results  in 

^  (A.12) 

Pi±lT'i*1  Wj-fl  - 

2  5  2  n 

A&-H  -b 


lVi=(G\ 
dt  \pj{ 


(A.13) 


for  internal  grid  points,  i=2,3,...N-l.  The  notations  p*± i,  ,  fw  and  r  —i  are 


defined  as 


Pi  ~  Pi±l 


(A.14) 


r,  r,  +  r,-±i 

1  i±i  =  — - — . 
—  2 

The  boundary  conditions  at  ip  =  0  is  described  by 


(A.15) 


(A.16) 


The  boundary  condition  at  ip  =  oo  is  described  by 

—Jf  =  £-^\pi±}. r^i(^v  -  ViV-i)]  -  -^[pN^T's-xiipN  -  93a/-i)],  (A.17) 
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and 


d-py 

di 


(A.18) 


This  results  in  2N(NSPCt1)  equations  for  NSPC  chemical  species  and  N  nodes. 
The  resulting  system  of  equations  is  still  computationally  stilt  due  to  the  extreme 
differences  in  gas  phase  reactions  rates.  A  special  algorithm  called  TRANSEQI  [2] 
was  used  to  make  the  equations  more  tractable. 

TRANSEQI  was  developed  by  Eraslan  and  Brown  [2]  to  reduce  stiffness  of  the 
differential  equations  describing  reactive  flow  problems.  In  addition,  core  memory  is 
reduced  and  a  nonstiff  ODE  solver  can  be  employed. 

Estimated  ionic  species  profiles  are  used  as  initial  conditions  in  the  algorithm. 
The  transport  coefficients  and  molar  production  rates  are  calculated  for  each  specie 
at  each  nodal  point  based  on  these  initial  conditions.  The  transport  equation  for 
the  first  specie  is  solved  holding  the  temperature  and  all  other  specie  concentrations 
constant.  At  the  next  time  step,  the  generation  equation  for  the  same  specie  is 
solved.  The  temperature  and  all  other  species  concentrations  are  again  held  con¬ 
stant  throughout  the  integration.  The  resulting  Aip  from  the  transport  and  chemical 
generation  equations  are  added  to  the  input  <p.  This  gives  the  new  ip  value  that  is 
used  in  subsequent  integrations.  The  transport  and  generation  equations  are  then 
solved  in  the  same  manner  for  each  of  the  remaining  species.  Property  and  mixture 
coefficients  are  then  calculated  using  the  new  values  of  specie  concentrations.  The 
evaluation  of  property  and  mixture  coefficients  can  be  updated  periodically  rather 
than  at  the  end  of  each  time  step,  which  helps  reduce  the  number  of  computations 
executed  in  the  program.  This  methodology  is  repeated  until  all  time  derivatives 
approach  zero.  Stiffness  is  virtually  eliminated  since  only  one  specie  is  treated  at  any 
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one  time  by  the  integrator.  Therefore,  Euler’s  scheme  and  Gear’s  method  are  used 

to  solve  Equations  (A.ll)  and  (A.13)  respectively. 

To  solve  Poisson’s  equation  (Eqn.  A.4)  it  was  rewritten  in  state  variable  form  by 

defining  the  electric  field  intensity  E  =  —  and  substituting  into  Equation  (A.4): 

dE  1  .  .  . .  , 

—  =  —  (n+  —  n_)e  (A.19) 

UX  6(3 


and 


dV 


(A.20) 

with  V(x=0)  and  V(x=L)  at  prescribed  boundary  conditions.  The  resulting  boundary 
value  problem  was  solved  using  a  shooting  method. 

The  predominate  ion,  cyclic  C3H3 ,  represents  99%  of  the  total  ion  current.  Since 
this  ion  is  treated  as  an  input  parameter  to  the  ionic  model,  Poisson’s  equation  (Eqn. 
A.4)  was  only  solved  at  the  end  of  the  first  iteration  in  TRANSEQI.  No  further 
changes  in  the  electric  field  are  expected  as  long  as  the  cyclic  C3H3  profile  remains 
fixed  [5]. 


B.  ESTIMATION  OF  INITIAL  CONDITIONS 


Good  initial  estimates  of  ion  concentration  profiles  can  result  in  substantial  de¬ 
creases  in  computation  time.  These  initial  profiles  were  obtained  by  neglecting  trans¬ 
port  effects  and  solving  the  resulting  plug  flow  model  described  by 

dYi  _  iOi 
dx  pu 

A  value  for  the  flame  velocity,  u,  can  be  taken  from  experimental  data  for  a  given  gas 
mixture  equivalence  ratio  and  pressure.  These  ionic  species  profiles  were  then  used 
as  initial  conditions  in  Equations  (A.ll)  and  (A.13). 
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III.  RESULTS 


The  results  from  the  plug-flow  model  were  compared  to  calculated  ion  profiles 
from  Frenklach.  When  using  Frenklach’s  neutral  species  profiles  the  plug-flow  model 
calculated  concentrations  of  major  ions  to  within  an  order  of  magnitude  of  the  values 
calculated  by  Frenklach’s  model.  Both  sets  of  calculated  profiles  exhibited  peaks  in 
concentrations  which  were  much  sharper  than  experimental.  The  calculated  plug-flow 
ion  concentrations  were  then  used  as  initial  conditions  in  the  diffusion  model  which 
eliminated  the  sharp  peaks  by  smoothing  and  broadening  the  profiles.  Since  the  peak 
concentration  did  not  change,  it  can  be  concluded  that  for  this  reaction  mechanism, 
the  plug  flow  model  can  be  used  to  give  good  quantitative  peak  concentration  values. 

The  use  of  ion-induced  electric  fields  to  calculate  ion  mobility  did  not  signif¬ 
icantly  alter  ion  profiles  compared  to  plug-flow  calculations.  Only  if  cyclic  C3H3 
becomes  a  variable  in  the  model,  rather  than  an  input  parameter,  can  ambipolar 
effects  become  apparent  as  ion  profiles  develop. 
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